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Abstract. We consider the GUT-like model with two scalar fields which has infinitesimal devi- 
ation from the conformal invariant fixed point at high energy region. In this case the dominating 
quantum effect is the conformal trace anomaly and the interaction between the anomaly-generated 
propagating conformal factor of the metric and the usual dimensional scalar field. This interaction 
leads to the renormalization group flow from the conformal point. In the super symmetric conformal 
invariant model such an effect produces a very weak violation of sypersymmetry at lower energies. 



1 Introduction 

43 



The formulation of consistent quantum field theory in curved space-time requires two new elements: 
the action of vacuum S v [g^] and also the nonminimal scalar-curvature interaction ^Rip 2 for 



each scalar field. The theory which lacks one of those elements is nonrenormalizable, because the 



corresponding counterterms appear anyway, and already at the one-loop level (see, for example, Q 
for the introduction). One can suppose that such a theory with an action S = S m i n + S non - m i n + S v 
should be valid until the effects of quantum gravity become important, that is for the energies 
well below the Planck scale. It is naturally to suppose that at the very high energies below the 
Planck scale the theory doesn't depend on the massive parameter and is conformal invariant. The 
condition of local conformal invariance is the special choice of the nonminimal parameter £ = ^ . 
The necessary action of vacuum has, in the conformal case, the restricted form: 

j4 



S v = J d 4 x^g~ ( ai C z + a 2 E + a 3 aR ) , (1) 

where \J —g C 2 = \/—g C^ va p C^ ua ^ is the conformal invariant square of the Weyl tensor and 
E = R 2 lua/3 — AR 2 lu J t-R 2 is the Gauss-Bonnet integrand. On quantum level the conformal invariance 
is spoiled by the trace anomaly which takes place in both vacuum j|] and matter-fields Q sectors 
of the theory. In the vacuum sector the trace anomaly arises (in the framework of dimensional 
regularization and the minimal subtraction scheme) because the counterterms, which must be added 

1 e-mail address: cognola@science.unitn.it 
2 e-mail address: shapiro@ibitipoca.fisica.ufjf.br 



1 



to the one-loop effective action, are local expressions in n-dimensional space-time. In general, the 
origin of the anomaly is the lack of regularization preserving both general covariance and local 
conformal invariance (see Q for the review of conformal anomaly and related issues). 

The interest to the conformal invariant theories in four dimensions (Ad) is partially inspired by 
the important role which they play in the two-dimensional (2d) theories which are traditionally 
considered as a toy models for the realistic higher dimensional theories. There are, however, a 
serious differences between 2d and Ad theories. First of all, in 2d the fields with a different spin 
contribute to the anomaly (which comes from the J d 2 x^—gR -type counterterm) with a different 
sign, that leads to the existence of the critical dimension for the sigma-models and thus in 2d the 
anomaly may be canceled. In Ad there are three possible types of the one-loop counterterms (|l|) 
and the fields with the different spin give contributions with the same sign to the renormalization 
constants of the vacuum parameters a\ and a 2 in (|l|). Therefore in Ad the critical dimension 
doesn't exist. The situation remains the same, even if the supersymmetry and quantum effects 
of gravity use to be incorporated Q. Another possibility which may be realized only in 2d is to 
start from the non-linear sigma-model in a background fields (defined as a geometric objects in 
target space) and thus provide the cancelation of anomaly in the noncritical dimension || 0. In 
Ad this scheme doesn't work either, because the nonlinear sigma-model is non-renormalizable and 
its possible higher derivative version H presumably contains unphysical massive ghosts. Thus, the 
quantum effects inevitably lead to the violation of conformal invariance through the trace anomaly 
with the consequent propagation of the conformal factor |], [lO|]p|. The study of the quantum theory 
of the conformal factor as a low-energy version of quantum gravity has been started in (See 
also 0). 

In the matter fields sector of the classically conformal invariant theory the one-loop divergences 
are also invariant but the finite one- loop contributions to the one- loop effective action are not. 
This can be, in particular, seen through the renormalization of the composite operators in the 
expression for the energy- momentum tensor Q . This renormalization leads to the violation of the 
conformal invariance in the higher-loop divergences |^]. Thus, in Ad, the classically conformal 
invariant theory suffers from two kinds of deceases: the conformal anomaly arises at one loop and 
produces the propagation of the conformal factor and also the nonconformal divergences take place 
at higher loops and break the renormalizability. Therefore the local conformal symmetry in Ad can 
not be exact, and can be realized only as an approximate high energy phenomena which may be 
called asymptotic conformal invariance. The asymptotic conformal invariance has been originally 



discovered in 16, 17] as a consequence of the conformal invariance of the one-loop divergences. 
This kind of the asymptotic conformal invariance doesn't follow from the asymptotic freedom and 
it puts some extra constraints on the multiplet composition and on the values of the coupling 



constants of the gauge model [18|. The shortcoming of this approach is that the higher order 
corrections to the /3-functions spoil the asymptotic conformal invariance and hence ^ = g is not 
the renormalization group fixed point beyond the one-loop level. Here we adopt another point 
of view on the asymptotic conformal invariance. Let us suppose that the asymptotically free or 
finite gauge theory in curved space-time, which is generated at Planck energy scale as an effective 
low-energy theory, is originally conformal invariant and thus the conformal invariant theory is the 
initial condition for the renormalization group flow in far UV rather than its fixed point. 

For the theories with a weak coupling the one-loop effects are always dominating and hence the 
leading quantum effect is the conformal (trace) anomaly 

< r/ >= k x C 2 + k 2 E + k 3 OR. 



1 The derivation of the anomaly-generated effective action of the conformal factor has been also performed in 
[ pd| for the case of the theory with torsion and in for the supersymmetric matter on the background of simple 
supergravity. 
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The values of k\,2 ,3 depend only on the number of fields of different spin in a GUT model The 
anomaly leads to the equation for the effective action 2g liu 5W/5g^ v = —^J—gT^. The solution of 
this equation is the Ad analog of the Polyakov action. It can be written in a local form []|] as: 

W[g^a] = S c [g^ u \ + j d 4 x^ { -aAa + a [k[C 2 + %(e- ~Dr\ } - A k' 3 R 2 \ . (2) 

The values of the coefficients k[ 2 3 are related to the ones of the trace anomaly in a well known 
way |J . The propagator of the conformal factor a is an inverse to the fourth derivative conformal 
invariant operator [19, ^, llCfl : 

A = n 2 + 2# fl/ V M V„ - \ru + \ ( V^ii) V M (3) 

and hence in the flat space limit a is massless field which can be essential at long distances. 
The solution (|2|) includes also an arbitrary conformal invariant functional S c . In four dimensional 
quantum field theory this functional can not be calculated exactly, and one can only establish its 
high energy asymptotic or derive some lower order terms in its expansion. However this functional 
doesn't depend on the conformal factor and thus it is irrelevant for our purposes 0. 

^From the physical point of view anomaly means the propagating conformal factor. Simul- 
taneously the two-loop effects produce the violation of the conformal constraint £ = g and as a 
result the conformal factor starts to interact with the usual scalar field. The investigation of the 
physical consequences of this interaction has been started in [24]. The interaction with the propa- 



gating conformal factor modifies the renormalization group equation for the scalar coupling / and 
£. These equations have minimal IR-stable fixed point £ = / = 0, and in the vicinity of this point 



one meets the first order phase transition, as a result the Einstein gravity is induced [24]. In [24] 
the theory with one scalar field was considered. Here we are going to perform the detailed study of 
the renormalization group equations in the framework of more complicated models with two scalar 
multiplets. Our purpose is to estimate the running of the couplings due to the interaction with 
the anomaly-generated conformal factor and, in particular, explore the renormalization group flow 
from the conformal fixed point in the models which are finite and supersymmetric in flat space 
time. 

Starting from the conformal initial point in far UV limit one can trace the renormalization 
group flows for the coupling constants and £ backward to the lower energies. In this way one 
can predict the value of £ at the lower energy scales. Thus we arrive at the necessity to study 
the renormalization group behavior in the theory which has only infinitesimal deviation from the 
conformal fixed point in far UV. To perform this in a consistent way one has to take into account 
the effect of conformal anomaly and the consequent quantum effects of the propagating conformal 



factor [24]. The detailed study of the renormalization group behavior is the purpose of the present 
paper. 

The paper is organized in a following way. In the next section we formulate the model of 
GUT coupled to the propagating conformal factor. In section 3 the derivation of the one-loop 
counterterms is performed. Section 4 is devoted to the renormalization group equations for two 
different models, and in the last section we draw our conclusions. 



2 The action of gravity induced by conformal anomaly has interest for the construction of quantum gravity Jl^, . 
It looks possible that the proper choice of the functional S c may revoke the discrepancy between the above effective 
action and direct calculations of W^^cr] in and with the results of the study of the two- and three-point 
functions El Etl. 
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2 Interaction of matter fields with conformal factor 



We suppose that the effects of conformal factor are relevant below the Planck scale where the 
nonminimal parameters £ become slightly different from 1/6 due to the higher loop effects. In 
curved space-time the transition to the low energies (or long distances) corresponds to conformal 
transformation in the induced gravity action jl], after that classical fields and induced gravity 
appear in a different conformal points ||24| (see also |13[). In order to take this into account one 
has to make a conformal transformation of the metric in (|2|) and then consider the unified theory. 
However it is more useful to perform the conformal transformation in the action of the matter 
fields. Such a transformation corresponds to some change of variables in the path integral for the 
unified theory, and doesn't modify the results of quantum calculations in the matter field sector, 
which we are interested in. 

The only source of conformal noninvariance in the action of the massless GUT model is the 
nonminimal term in the scalar sector. As far as the values of the parameters £ are not equal to | the 
conformal factor starts to interact with the scalar fields. Consider the general SU(N) model with 
two types of scalar fields: real ones $ a in the adjoint representation and complex ones ip l in vector 
representation of the gauge group. In curved space-time one has to introduce two nonminimal 
parameters £1 and £2 ~ one for each type of the scalar fields. This model has been investigated 



in [25] without taking into account the back reaction of vacuum, and one can find more complete 
information including the full set of the /^-functions in this paper. Below we write down only that 
parts of the /3-functions which we shall actually use. 

Introducing the scale parameter a we arrive at the following action for the conformal factor 
coupled to scalar fields: 



S sc = W[ gfiu ,a}+ / (fx 



~ (1 - 6^)$^ + (1 - 6£ 2 )^V 



(a 2 (Vcr) 2 + a Oct) 



+~ sr {v^T {v u ^) a + sr {v^ + )i (Putf + \ ixR^ a + £ 2 i^V - v($ a , ^)| , (4) 

where Wfg^cr] has been defined in (||) and we use the notation (Vex) 2 = g^ v d^odyO . In all the 
expressions T> are the derivatives of the matter fields which are covariant with respect to both 
gauge and gravitational fields. The notations for the SU(N) group including the relations between 

symmetric D ra b and antisymmetric f abc structure constants, generators (jj^ , traces etc. can be 
found in pi (see also second reference in p6fl ) . The flat-space part of the potential has the form 



V(<S> a , tf) =±h ( W) 2 + i h ($ a D rab $ b ) 2 + \h (W) (^V) 



(5) 



The action of scalar fields (f|) must be supplemented by the action of spinors and gauge fields which 
are part of the GUT model. The corresponding Lagrangian has the form [25] 



1 m m / \ a \3 

£ = — A {G%? + i £ * a ik) {fvf - hif*>&) *fo + • £ i,\ k) [ycj? -h 2 ( — ) ^ m + 

4 k=l k=l V 1 J i 



s=l 



k=l 



X a \ J 



(6) 
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We are interested in the quantum theory of matter fields and conformal factor a on the back- 
ground of the classical metric. The interaction between matter fields and conformal factor arises as 
a result of the conformal transformation of the metric g^ u — > g' = g^ v exp(2acr) and the matter 
field w — > w' = w exp(d w aa), where d w is the conformal weight of the field w. One can see that 
the massless spin-^ and spin-1 fields decouple from the conformal factor. The only kind of fields 
which takes part in such an interaction are scalars, for which the interaction with conformal factor 
appears as a result of the violation of the condition £ = 1/6 at low energies. The contributions of 
other matter fields to the /3-functions and effective potential of the scalar fields do not depend on 
the conformal factor, they are given by the usual expressions |25|, [l]]. 



3 Calculation of the one-loop divergences 

The calculation of the one-loop divergences of the theory (|j) is not a trivial problem because the 
classical action contains second derivative terms as well as fourth derivative ones. Here we shall use 
the approach of Ref. |26| , where the one- loop divergences have been calculated for higher derivative 
quantum gravity coupled to the matter fields. Following [26|, we shall use the background field 



method and generalized Schwinger-DeWitt technique of 27]. Thus we start from the separation of 



fields into background a, 
where 



and quantum r, r\ ones by changing (a, <3? a , <p % , c/\ 



<<P ,<Pj 



a 



a + T , $ a/ = $ a + irf , 



<P+i)C, <Pi 



<Pi +iXi 



(7) 



and the imaginary units are introduced for convenience. The one-loop effective action is defined as 



Tr In H, 



(8) 



where H is the bilinear (with respect to the quantum fields r, rj, x, X + ) form of the classical action 
(Eh. After some algebra we get the following self-adjoint bilinear form 



H 



( H TT 
H, 
H 



?/r 



H 
H 
H 
H 



TT) 

nv 



X'l 



H x+x H x+x+ 
H xx H 



x\ n 



0) 



with the following derivative structure of the operator H: 



H 



(10) 



The bilinear form of the action is matrix differential operator with fourth derivatives in the H TT 
sector and with second derivatives in the sector of usual scalar fields and in the mixed pieces. This 
structure of the bilinear form is similar to the one which is known from the theory of multiscalar 



GUT coupled to higher derivative gravity 26]. This is indeed natural, because what we are doing 
now is nothing but the study of induced quantum gravity (Q) unified with the same GUT model. 
This analogy facilitates the calculations considerably, because of the following three reasons: 
i) The general expression for the divergent part of | Tr In H ( |i~0| ) is known from |26[] 



n— 4 



2 / d n xyf=g tr 



ippQi + Ya v>xvv ' 



[IV 



+ 



1 

48 



ra 2 



1 

6 



V^R UU + R - U 



12 



5 



1 



1 



+ __ IRTR^ --^j- — R 2 | + I Tr In {in + £ A V A + 1)}^ + (surface terms) . (11) 



The above formula includes the standard contribution of the second order operator 



-Tr In {!□ +E x V x + d] 



div 



(see, for example, the contribution of the minimal fourth order operator first calculated in 
]i~0[ | and also the contributions from the mixed sector derived in Ref. [26| for the case of higher 
derivative gravity coupled to matter fields. The formula ( |i~l"l) shows, in particular, that we do not 
need explicit expression for N x . 

ii) Calculating the divergences with the use of the formula (|TT|) one can learn the form of the 



beta- functions. However, since the formula (11) is the same as for the higher derivative gravity, it 
is easy to see that the general structure of the renormalization group equations in our theory is 
also the same as the one established in pq| . In particular, the quantum gravity corrections to the 
beta-functions of the nonminimal parameters are completely universal, they do not depend on 
the model and have the same form as for the one-scalar model [24]. Therefore all that we need 
is the form of the quantum gravity corrections to the beta- functions of the constants f%. These 
corrections can be calculated on flat background and hence in what follows we put g^ v = r]^ u and 
switch off all the curvature dependent terms. 



iii) The structure of the expression (11) is direct generalization of the one we have already studied 
in (24]] for the 1-scalar model. One can easily check that here, just as in the 1-scalar model, the 
contributions of tv PiQi, tr (P1Q2 + P2Q1), tr (P1Q3 + P3Q1), and tr P2Q2 cancel, and therefore 
we need only the expressions for V 1 *", U, Pi, Q%, E x , D. Disregarding the curvature dependent 
terms one can obtain the following expressions for V, U, Qi, Pi, D ■ 



Pi = —ia 



yv = -a 2 (Ci <S> a $ a + C2 vttf) 9^ , U = 0, 

Qi = -ia (Ci$ 6 (2^ C 2 ^ + ) , 



V J 



-CiZ5 ab + y ab Vf y a i 

D = H=\ Vj -(2Z6} + v? V** I , (12) 



3 1 
'ij 



Vib Vij ( 2 zsl + yl 



where Z = a 2 (Vcr) 2 + a(D<r) and 



ir \ 1 



I = A$2 5 a6 + hq>a ^b + /2 (\ Drah D rcA + D rac D rbd ) <S> c <S> d + f 3 5 ab <p 2 + UD rab iff ( y ) ^ , 



j 



V aj = / 3 *V/ + UDrac^t ( y) > Vib = /s* V + UDrbc$ C <pt (y 

Vj = / 3 <I>V + hD rac ^ ( £Y /, vi = / 3 <&v + UD rac ^ (£X 



2 



k 



1 1 / A r \ 1 

ri _ f_rtsCfi.cn 1 ^ r fi,cr> I * \ , f , 2;i 



VJ = -h^ c 8) + -h<S> c D rcd § a + h^5) + f B <pt<p 

V"-' V ij = f 5 <p+cp+. 
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Substituting these expressions into ([Tl]) we arrive at the explicit form of the corrections to the 
matter sector of T^iv from the quantum conformal factor^]. 



r (D = f^l 

div g. 



2 J ^ x ^{ia 4 (l-Ci) 2 (W) 2 + 2^ 



+2a 4 CiC 2 (l - Ci - C2 + CiC 2 )$ V*> + lf2a\l(^ a D abc ^ b f + / 3 a 2 (( + 4CiC 2 + (2 &<P + <P 



-2 , a/- /■ 1 /~2 



,2, 



+ha 2 [Ci 2 + 4CiC2 + C 2 2 ] {<5> a D abc <5> h )vt (y)V' + 6/5a 2 C 2 2 (^ + ^ 2 | , (13) 

where e = (47r) 2 (n — 4) is the parameter of dimensional regularization, and we disregarded all 
surface and matter-independent terms. 

One can see that there are no any divergences which lead to renormalization of scalar fields 
<£, ip. Thus the only modifications due to the contributions of the conformal factor are in the 
renormalization of couplings /12 .. 5 and parameters £12, and the renormalization group equations 
include only (3— functions, but not 7-functions. The contribution of the conformal factor to the 
effective potential of the scalar fields depends on the /3-functions only. 

4 Renormalization group equations 

Now we are in a position to perform the renormalization group study of the back-reaction of vacuum 
to the matter fields. For our purposes it is convenient to use the formulation of the renormalization 
group in curved space-time, given in [[17], The renormalization group equation for the effective 
action in curved space-time has the form 

{^djj, + /3p ^p ~ I dAy f~y iw sJ( x ) } T \ w > p > 9 ^>^ = °> ( 14 ) 

where w is the full set of the quantum fields (gauge, spinor, Higgs and conformal factor a) and 
p - complete set of parameters including £'s. The solution corresponding to the desirable scaling 
behavior has the form 

^[e^g^u, w, p, fi] = Tlg^, w(t), p(t), fi\ , (15) 

where fi is the dimensional parameter of renormalization. Effective fields and coupling constants 
obey the equations 

(4vr) 2 ^^ = ( lw + d w )ct>, w(0) = w, (16) 

at 



(^) 2j ^ l = P P +pd P , P(0)=p, (17) 

where 7 and functions are defined as usual. According to the results of the previous section the 
renormalization of the fields and all couplings except fi and £j are not modified by the contributions 
of the quantum field a. Thus the /^-functions for the gauge and Yukawa coupling constants are 
just the ones derived in [^]. On the other hand, these contributions to the j3f i and (3^ all have 

3 The renormalization in the vacuum sector was discussed in ^] for the case of the one-scalar model. Its form 
doesn't depend on the number of scalar fields and that is why we do not discuss it here. 
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universal form and do not depend on the gauge group of the theory. Thus we find for our theory 
(from this moment we will be using the variables Q = 1 — 6£j for convenience): 

fa = $ + A ^ . ^ = 4 0) + ' = $t > & = $ 0) < ( 18 ) 

where (3^ is the /3-function for the effective parameter p in curved space-time without back reaction 
of vacuum (or other form of quantum gravity) and A(5 P are the quantum gravity corrections. In 
our case, contrary to the high derivative gravity |^] those A/? p 's are nonzero only for the scalar 
and the nonminimal parameters. 

According to recent communications (see, for example, [pg|] ) the values of £j are very important 
at the energies between the Planck scale Mpi = 10 19 CreV and the unification scale Mx = 10 14 GeV, 
because there is a hope to meet natural inflation for this rate of character energies. Thus our 
purpose is to evaluate the running of £j (or Q) backward from the high energy Planck scale. As 
a lower end of the energy interval one can take the unification point Mx = 10 14 ~GeV but, since it 
doesn't lead to the serious changes in the calculations we shall take, as a lower limit, the Fermi 
scale Mp = 100 GeV. 



4.1 SU(N) model: running away from the finite theory 

As an examples of the effect of the quantum conformal factor we shall consider two toy models 



with two scalar fields: one which possesses the one-loop finiteness without supersymmetry [25], and 
another one which admits N = 2 supersymmetry but is not finite. 

Let us start with the generic model (Q) . The counterterms ( |l3|) lead to the following expressions 
for Af3 p : 



(4vr) 2 A(3 fl = 12a 2 h (( + 4a 4 d 2 (Ci - I) 2 , /i(0) = fi ; (19) 



(4vr)^A/3 /2 = 12a 2 /2Ci 2 



/ 2 (0) = / 2 ; (20) 



(4tt) 2 A/3/3 = 2a 2 h (Ci + 4CiC 2 + Cl) + 4a 4 CiC 2 (1 - Cl - C2 + C1C2) , / 3 (0) = h \ (21) 



(4vr) 2 A/3 /4 = 2a 2 / 4 (C? + 4CiC 2 + Cl) , 

(4vr) 2 A/3 /5 = 12a 2 h Cl + 4a 4 Cl(C 2 - I) 2 , 
and 



/ 4 (0) = / 4 ; (22) 
/ 5 (0) = h (23) 



(4^) 2 A/3 C . = 2a 2 C 2 (Ci - 1) , 0(0) = d • (24) 

Let us now consider the asymptotic behavior of the effective couplings fi(t),Q(t). As it was 
already mentioned above, all the /3-functions (corresponding to gauge, Yukawa, scalar couplings 
and to the nonminimal parameters £) vanish in the conformal fixed point 

h k {t) = h%, fi(t) = f*, = 0, (25) 

where /*, h\ are the values corresponding to the fixed point in flat space-time. The corrections 
Ap indicate that there is also a second "minimal" fixed point with the same (p5|) solutions for h^{€) 
and fi(t) but with Q = 1 (this corresponds to the £j = in (Q), that is why we call this fixed 
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point minimal). The behavior of the effective charges and effective potential in the vicinity of the 



minimal fixed point has been studied in [24 for the one-scalar model. Now we are interested in the 
behavior of the effective charges close to the conformal fixed point. 

As far as the renormalization group equations for g(t), hk(t), fi(t), Cj(t) are very cumbersome, 
it is reasonable to take particular models in which the study of these equations performs easier. 
As a first example, let us take a particular model of (||) with N = 8 and one scalar multiplet in 
the adjoint representation. When the number of spinor multiplets is m = 84 the theory is one-loop 
finite in the flat space-time (see Ref. [^] for a detailed discussion of such a model without taking 
the quantum effects of the conformal factor into account). To perform the numerical analysis of the 
renormalization group equations one has to impose the initial conditions at far UV. We suppose that 
the running from the conformal fixed point happens because of the two-loop contributions which 
violate the £ = | condition and thus switch on the interaction with the conformal factor. Thus, to 
impose the initial conditions one has to evaluate the two-loop effects. The two-loop contributions 

4 

are, generally speaking, proportional to multiplied by the combinatorial factor. Let us take 
the typical value of the parameter £ as 10 -3 , while the coupling constant is g 2 = 0.1. Then the 
numerical analysis enables one to trace the behavior of the effective couplings C(t), f 12(f) for 
> t = ln(fi/Mpi) > —39 that corresponds to the running between fi = M Pi = 10 19 GeV and 
fj, = Mp = 10 2 GeV. 

We have performed the numerical analysis of the full system of the renormalization group 
equations for the effective gauge, Yukawa, scalar couplings and £(£) with taking into account the 
corrections (|H]) - fl24]). As an initial point of the renormalization group flow in far UV we choose 



the values (see [25] for the details) 

h 2 = 0; ftm 0.43621 , / 2 * ps 0.30277 

that provides the one-loop finiteness in the theory without back reaction of vacuum. The results 
for the contribution to the couplings due to back-reaction, that is the differences Dfj and DQ 
between the couplings in the with and without the back-reaction of vacuum, are plotted in Fig. |l|. 
It is easy to see that the plots for the deviations look like a plane lines. The reason for this is 
that the effect of the back reaction is very week and in such a "small" interval on the logarithmic 
scale the non-linear function looks like a linear. Also the numerical values presented at the Fig. [l] 
show that the effect is very small. For instance, the deviation of /1 2 from the finite fixed point is 
about six orders smaller than the value of the coupling itself in this fixed point. One has to notice 
that since the effect is almost linear, it is quite easy to construct the same plots for other initial 
deviations DQ\, 

The similar numerical analysis has been performed for the more interesting SU (5) model (||) - 
(rc3|) with two kinds of scalar fields and with 



m = l, n=15; hl 23 (t) = h* 123 , ft.,.,S) = /i*..,5 ; 



where 



^1,2,3 = ( L421 , 1-681, 2.361) ; /i*2,3,4,5 = (0.659, 1.293, 0.324, 1.677, 1.039). 

The results of the numerical solution are presented at the Fig. || .The parameter DQ\(Q) here is 
taken 0.01 as in the previous example, and we kept the constraint DQ\ = D( 2 for the sake of 
simplicity. Qualitatively the results are the same as in the one-scalar case - the effect exists but it 
is tiny. 
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Figure 1: Contribution to the running coupling constants due to back-reaction. SU [8] model with one 
scalar field and 84 spinor multiplets. 




4.2 Model with the rigid breaking of supersymmetry 

Let us now consider the SU(2) gauge model with an action |2^, |16|, |l[: 

+i ip a 7 M V° b ip b - hii> a F abc (f b 1p c -i h 2 ^al5 Fabc A b ip c + 



2\2 



h(A 



2\2 



h<p 2 A 2 



(<pAf 



(26) 



where (p and A are scalar and pseudoscalar fields in the adjont representation of the SU (2) gauge 
group. Other notations are obvious. In flat space, for the particular values of the couplings 



h = f2 = h = 0, h = h\ = hl = g 2 



(27) 



the model (|26| ) possesses the N = 2 supersymmetry. For other values of the couplings one meets 
the renormalizable theory which has also second non-supersymmetric renormalization group fixed 
point |29| . In curved space time this model is also renormalizable (if only the nonminimal terms 
and the action of vacuum are introduced), and admits the asymptotic conformal invariance for (f2^) 
and £i + & = 1/6. 
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Taking into account our previous results we arrive at the renormalization group equations for 
the above model with the back reaction of vacuum. The gauge coupling behaves as |2£ 



2, * - 1 



and supersymmetry fixes Yukawa constants to be h\ = h% = g 2 . In flat space-time those relations 
hold under renormalization p9fl . The same happens when we put our theory in curved space- 
time and take the back reaction of vacuum into account, because the quantum conformal factor 
contributes to the /^-functions of the scalar couplings and nonminimal parameters. We are interested 
in the renormalization group flow from the supersymmetric and conformal invariant UV fixed point, 



that is why we have to impose the constraints (27) on the initial data. As far as the back reaction 



doesn't concern the behavior of the Yukawa couplings, the relations (|27]) for the Yukawa couplings 
hold for all scales and they can be indeed used in the RG equations directly, while 71,2,3.4 should be 
regarded as arbitrary quantities. After we use the relations between gauge and Yukawa coupling 
constants, the RG equations for /1, 2,3,4 (t) become 

(4vr) 2 ^ = lift + 3/! + 8/3/4 + 8/| - 8f l9 2 - 8g* + 12a 2 /iCi 2 + 4a 4 C 2 (Ci - I) 2 , 

(4n) 2 ^ = + 3/| + 8/3/4 + 8/| - 8f 2 g 2 - 8g 4 + 12a 2 / 2 C 2 2 + 4q 4 C 2 2 (C 2 - I) 2 , 

(4tt) 2 ^ = 4/| + (5/3 + 4/ 4 )(/i + / 2 ) + 8/| - 8/ 35 2 - 8 5 4 + 2a 2 / 3 (d 2 + 4CiC 2 + C 2 2 )+ 
+ 8a 4 CiC2(l-Ci-C2 + CiC 2 ), 

(4vr) 2 ^ = 6/| - 8/ 49 2 - 6 5 4 + 8/3/4 + 2/ 4 (/i + h) + 2a 2 / 4 (d 2 + 4CiC 2 + C 2 2 ) , (28) 
while for £1,2 (/) we meet the equations: 

(4vr) 2 ^ = Ci (5/i - 4 5 2 ) + C 2 (3/ 3 + 4/4) + 2a 2 C 2 (Ci - 1) , 



(4vr) 2 = (2 (5/2 - 4 5 2 ) + Ci(3/ 3 + 4/4) + 2a 2 C 2 2 (C2 - 1) . (29) 

In order to study the renormalization group behaviour of the effective couplings f(t), £(t) we 
have to choose the initial value of g 2 and the initial variation DC,\ = —D^ 7^ 0. As in the previous 
case we take in the far UV g 2 = 0.1 and Ci = 0.001 . The results of the numerical analysis, 
concerning the contributions due to back-reaction, are plotted in Figs. |3| for the case £i = 0.001. 

One can see that in this case the effect of quantum conformal factor is essentially the same as in 
the previous case of the SU(8) and SU(5) models. Because of those effects one meets the nonzero 
/i,2,3 couplings with the first two having the values about 10~ 8 at the Mx unification scale and 
about 10~ 6 at the Mp Fermi scale, while they remain identically zero without the back reaction of 
the conformal factor. 

In order to give some speculation on the strongly coupled theories we also investigated the case 
with initial values of g 2 = 0.5 and £ = 0.1. The results are plotted at Fig. |||. 
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Figure 3: Contribution of back-reaction to the running couplings fj(t) and Q(t). First example with 

Ci(o) = -C 2 (o) = o.ooi. 




Figure 4: Contribution of back-reaction to the running couplings fj(t) and (j(t). g 2 — 0.5 and Ci(0) = 

-C 2 (o) = o.i. 

5 Conclusion 

The arguments were presented that in four-dimensional quantum field theory in curved space-time 
the local conformal symmetry can not be exact, and that it can be only approximate symmetry 
which serves as the initial point for the renormalization group equations in the far UV limit. Starting 
from the conformal invariant theory one meets the trace anomaly and consequent propagation of 
the conformal factor of the metric. We have considered the interaction between the quantized 
conformal factor and the matter fields in a theory with approximate conformal symmetry, in a 
region close to the scale of asymptotic freedom and asymptotic conformal invariance. 

The contributions of the conformal factor to the /3-functions of the scalar coupling constants 
have been calculated and it was shown that these contributions drive the theory out from the 
conformal point and also modify a values of the scalar coupling constants. If the starting model in 
far UV possesses supersymmetry and (or) finiteness, the interaction with the conformal factor leads 
to the violation of those properties at lower energies. For the initial violations about DQ = 0.001 
the numerical values of the deviations of the scalar couplings from the symmetric state range from 
10~ 8 to 10~ 6 at the Fermi scale and about one order less at the Unification scale. Indeed the 
effect is very weak and therefore all the dependencies are very close to be linear. In particular the 
above deviations may be considered as proportional to the squares of the initial violations (DQ) 2 
of the conformal invariance. At higher loops the supersymmetry breaking in the scalar sector 
should violate supersymmetry in the Yukawa and gauge interactions. As a result the relations 
between the masses of the particles at low energies may differ from the ones which are expected 
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Figure 5: Contribution of back-reaction to the running couplings fj(t) and Cj{t)- 9 2 — 0-5 and Ci(0) = 
— C2(0) = 0.5. This example was included only to illustrate the nonrealistic order of magnitude for g 2 and 
initial deviations DQ for which the nonlinear behavior shows up. 

in supersymmetric GUT's and the supersymmetry is not seen at low energies. However, since the 
numerical effect of such a violation will be extremely small, it is very difficult to see how one can 
observe it. One can say that if the local conformal invariance exists as a high energy symmetry, it 
holds as a very good approximation at lower energies. 
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Back reaction of vacuum and the renormalization group How 
from the conformal fixed point 
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Abstract. We consider the GUT-like model with two scalar fields which has infinitesimal deviation 
from the conformal invariant fixed point at high energy region. In this case the dominating quantum 
effect is the conformal trace anomaly and the interaction between the anomaly-generated propagating 
conformal factor of the metric and the usual dimensional scalar field. This interaction leads to the 
renormalization group flow from the conformal point. In the supersymmetric conformal invariant 
model such an effect produces a very weak violation of sypersymmetry at lower energies. 

1 Introduction 

The formulation of consistent quantum field theory in curved space-time requires two new elements: 
the action of vacuum Sylg^] and also the nonminimal scalar-curvature interaction ^Rtp 2 for 
each scalar field. The theory which lacks one of those elements is nonrenormalizable, because the 
corresponding counterterms appear anyway, and already at the one-loop level (see, for example, [1] 
for the introduction) . One can suppose that such a theory with an action S = S m ,„ + S non - m i n + S v 
should be valid until the effects of quantum gravity become important, that is for the energies well 
below the Planck scale. It is naturally to suppose that at the very high energies below the Planck 
scale the theory doesn't depend on the massive parameter and is conformal invariant. The condition 
of local conformal invariance is the special choice of the nonminimal parameter £ = ^ . The necessary 
action of vacuum has, in the conformal case, the restricted form: 

S v = J d 4 Xv ^( ai C 2 + a 2 E + a 3 OR) , (1) 

where yj—g C 2 = yj—g C^ va fj C^ va ^ is the conformal invariant square of the Weyl tensor and 
E = R 2 liya p — 4i? 2 (J/ + R 2 is the Gauss-Bonnet integrand. On quantum level the conformal invariance 
is spoiled by the trace anomaly which takes place in both vacuum [2] and matter-fields [3] sectors 
of the theory. In the vacuum sector the trace anomaly arises (in the framework of dimensional 
regularization and the minimal subtraction scheme) because the counterterms, which must be added 
to the one-loop effective action, are local expressions in ra-dimensional space-time. In general, the 
origin of the anomaly is the lack of regularization preserving both general covariance and local 
conformal invariance (see [4] for the review of conformal anomaly and related issues) . 

1 e-mail address: cognola@science.unitn.it 
2 e-mail address: sliapiro@ibitipoca.fisica.ufjf.br 
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The interest to the conformal invariant theories in four dimensions (4c?) is partially inspired 
by the important role which they play in the two-dimensional (2c?) theories which are traditionally 
considered as a toy models for the realistic higher dimensional theories. There are, however, a serious 
differences between 2c? and 4c? theories. First of all, in 2c? the fields with a different spin contribute 
to the anomaly (which comes from the J (PxyJ—gR-type counterterm) with a different sign, that 
leads to the existence of the critical dimension for the sigma-models and thus in 2c? the anomaly 
may be canceled. In 4c? there are three possible types of the one-loop counterterms (1) and the 
fields with the different spin give contributions with the same sign to the renormalization constants 
of the vacuum parameters a\ and a 2 in (1). Therefore in 4c? the critical dimension doesn't exist. 
The situation remains the same, even if the supersymmetry and quantum effects of gravity use to be 
incorporated [5]. Another possibility which may be realized only in 2c? is to start from the non-linear 
sigma-model in a background fields (defined as a geometric objects in target space) and thus provide 
the cancelation of anomaly in the noncritical dimension [6, 7]. In 4c? this scheme doesn't work either, 
because the nonlinear sigma-model is non-renormalizable and its possible higher derivative version 
[8] presumably contains unphysical massive ghosts. Thus, the quantum effects inevitably lead to the 
violation of conformal invariance through the trace anomaly with the consequent propagation of the 
conformal factor [9, 10] 1 . The study of the quantum theory of the conformal factor as a low-energy 
version of quantum gravity has been started in [13] (See also [14]). 

In the matter fields sector of the classically conformal invariant theory the one-loop divergences are 
also invariant but the finite one-loop contributions to the one-loop effective action are not. This can 
be, in particular, seen through the renormalization of the composite operators in the expression for the 
energy-momentum tensor [3]. This renormalization leads to the violation of the conformal invariance 
in the higher-loop divergences [15]. Thus, in 4c?, the classically conformal invariant theory suffers 
from two kinds of deceases: the conformal anomaly arises at one loop and produces the propagation of 
the conformal factor and also the nonconformal divergences take place at higher loops and break the 
renormalizability. Therefore the local conformal symmetry in 4c? can not be exact, and can be realized 
only as an approximate high energy phenomena which may be called asymptotic conformal invariance. 
The asymptotic conformal invariance has been originally discovered in [16, 17] as a consequence of the 
conformal invariance of the one-loop divergences. This kind of the asymptotic conformal invariance 
doesn't follow from the asymptotic freedom and it puts some extra constraints on the multiplet 
composition and on the values of the coupling constants of the gauge model [18]. The shortcoming of 
this approach is that the higher order corrections to the /3-functions spoil the asymptotic conformal 
invariance and hence £ = g is not the renormalization group fixed point beyond the one-loop level. 
Here we adopt another point of view on the asymptotic conformal invariance. Let us suppose that the 
asymptotically free or finite gauge theory in curved space-time, which is generated at Planck energy 
scale as an effective low-energy theory, is originally conformal invariant and thus the conformal 
invariant theory is the initial condition for the renormalization group flow in far UV rather than its 
fixed point. 

For the theories with a weak coupling the one-loop effects are always dominating and hence the 
leading quantum effect is the conformal (trace) anomaly 

< T/ >= k x C 2 + k 2 E + k 3 nR. 

The values of k\, 2 ,3 depend only on the number of fields of different spin in a GUT model [2]. The 
anomaly leads to the equation for the effective action 2g flu SW/Sg^ = —yj—g Tjf. The solution of 
this equation is the 4c? analog of the Polyakov action. It can be written in a local form [9] as: 

W[g„, a] = S c [g^} + J j ^aAa + a [k[C 2 + k' 2 (e - jnR^j ] - ± k' 3 R 2 } . (2) 

1 The derivation of the anomaly-generated effective action of the conformal factor has been also performed in [11] for 
the case of the theory with torsion and in [12] for the supersymmetric matter on the background of simple supergravity. 
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The values of the coefficients k[ 2 3 are related to the ones of the trace anomaly in a well known 
way [9]. The propagator of the conformal factor a is an inverse to the fourth derivative conformal 
invariant operator [f9, 9, 10]: 

A = n 2 + 2i?^V A1 V i ,-^i?n + ^(V^)V A1 (3) 

and hence in the flat space limit a is massless field which can be essential at long distances. The solu- 
tion (2) includes also an arbitrary conformal invariant functional S c . In four dimensional quantum 
field theory this functional can not be calculated exactly, and one can only establish its high en- 
ergy asymptotic or derive some lower order terms in its expansion. However this functional doesn't 
depend on the conformal factor and thus it is irrelevant for our purposes 2 . 

From the physical point of view anomaly means the propagating conformal factor. Simultan- 
eously the two-loop effects produce the violation of the conformal constraint £ = g and as a result 
the conformal factor starts to interact with the usual scalar field. The investigation of the physical 
consequences of this interaction has been started in [24]. The interaction with the propagating con- 
formal factor modifies the renormalization group equation for the scalar coupling / and £. These 
equations have minimal IR-stable fixed point £ = / = 0, and in the vicinity of this point one meets 
the first order phase transition, as a result the Einstein gravity is induced [24]. In [24] the theory with 
one scalar field was considered. Here we are going to perform the detailed study of the renormaliza- 
tion group equations in the framework of more complicated models with two scalar multiplets. Our 
purpose is to estimate the running of the couplings due to the interaction with the anomaly-generated 
conformal factor and, in particular, explore the renormalization group flow from the conformal fixed 
point in the models which are finite and supersymmetric in flat space time. 

Starting from the conformal initial point in far UV limit one can trace the renormalization group 
flows for the coupling constants and £ backward to the lower energies. In this way one can predict 
the value of £ at the lower energy scales. Thus we arrive at the necessity to study the renormalization 
group behavior in the theory which has only infinitesimal deviation from the conformal fixed point 
in far UV. To perform this in a consistent way one has to take into account the effect of conformal 
anomaly and the consequent quantum effects of the propagating conformal factor [24]. The detailed 
study of the renormalization group behavior is the purpose of the present paper. 

The paper is organized in a following way. In the next section we formulate the model of GUT 
coupled to the propagating conformal factor. In section 3 the derivation of the one-loop counterterms 
is performed. Section 4 is devoted to the renormalization group equations for two different models, 
and in the last section we draw our conclusions. 

2 Interaction of matter fields with conformal factor 

We suppose that the effects of conformal factor are relevant below the Planck scale where the non- 
minimal parameters £ become slightly different from 1/6 due to the higher loop effects. In curved 
space-time the transition to the low energies (or long distances) corresponds to conformal transform- 
ation in the induced gravity action [1], after that classical fields and induced gravity appear in a 
different conformal points [24] (see also [13]). In order to take this into account one has to make 
a conformal transformation of the metric in (2) and then consider the unified theory. However it 
is more useful to perform the conformal transformation in the action of the matter fields. Such a 
transformation corresponds to some change of variables in the path integral for the unified theory, and 

2 The action of gravity induced by conformal anomaly has interest for the construction of quantum gravity [14, 20]. 
It looks possible that the proper choice of the functional S c may revoke the discrepancy between the above effective 
action and direct calculations of W\_g^ V) a] in [21] and with the results of the study of the two- and three-point functions 
[22, 23]. 
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doesn't modify the results of quantum calculations in the matter field sector, which we are interested 
in. 

The only source of conformal noninvariance in the action of the massless GUT model is the 
nonminimal term in the scalar sector. As far as the values of the parameters £ are not equal to ^ 
the conformal factor starts to interact with the scalar fields. Consider the general SU(N) model 
with two types of scalar fields: real ones $ a in the adjoint representation and complex ones ip l in 
vector representation of the gauge group. In curved space-time one has to introduce two nonminimal 
parameters £1 and £2 ~~ one for each type of the scalar fields. This model has been investigated 
in [25] without taking into account the back reaction of vacuum, and one can find more complete 
information including the full set of the /3-functions in this paper. Below we write down only that 
parts of the /3-functions which we shall actually use. 

Introducing the scale parameter a we arrive at the following action for the conformal factor 
coupled to scalar fields: 



S, 



cr] + J d 4 2 
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\(\- 66)$ a $ a + (1 - 



a (Vo-) 2 + ana) 



+± <r (iw (v^y + <r {v^y (v^y + ±& i? w +6^V - n$w)} , (4) 

where W\g^ v ^a\ has been defined in (2) and we use the notation (V<r) 2 = g^ v d^ad v a . In all the 
expressions V are the derivatives of the matter fields which are covariant with respect to both 
gauge and gravitational fields. The notations for the SU(N) group including the relations between 

symmetric D ra b and antisymmetric f a i c structure constants, generators (r^j , traces etc. can be 
found in [I] (see also second reference in [26]). The flat-space part of the potential has the form 

tf) = i f x (W) 2 + i f 2 {<5> a D rah <5> h ) 2 + \ h (^V) 



+\u{^D rab ^)^t (£|V + ^/ 5 (^V) 2 . (5) 

The action of scalar fields (4) must be supplemented by the action of spinors and gauge fields which 
are part of the GUT model. The corresponding Lagrangian has the form [25] 

k=i k=i v j 1 

s=l k=l J 

We are interested in the quantum theory of matter fields and conformal factor a on the background 
of the classical metric. The interaction between matter fields and conformal factor arises as a result 
of the conformal transformation of the metric g^ v — > g' = g^ v exp(2a<r) and the matter field 
w — > w' = w exp(d w aa), where d w is the conformal weight of the field w. One can see that the 
massless spin-^ and spin-I fields decouple from the conformal factor. The only kind of fields which 
takes part in such an interaction are scalars, for which the interaction with conformal factor appears 
as a result of the violation of the condition £ = 1/6 at low energies. The contributions of other matter 
fields to the /3-functions and effective potential of the scalar fields do not depend on the conformal 
factor, they are given by the usual expressions [25, I]. 
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3 Calculation of the one-loop divergences 



The calculation of the one-loop divergences of the theory (4) is not a trivial problem because the 
classical action contains second derivative terms as well as fourth derivative ones. Here we shall use 
the approach of Ref. [26], where the one-loop divergences have been calculated for higher derivative 
quantum gravity coupled to the matter fields. Following [26], we shall use the background field 
method and generalized Schwinger-DeWitt technique of [27]. Thus we start from the separation of 
fields into background a, 4> and quantum r, rj ones by changing (a, $ a , p l , (pf) — > (a 1 , $ a/ , <p 1 ' , cpf ) 
where 



a' = a + T, $ a/ = $ a + i?7 a , p 1 ' = (p l + i X \ <pf = <pf + ixf (?) 
and the imaginary units are introduced for convenience. The one-loop effective action is defined as 

r = ^Trlnif, (8) 

where H is the bilinear (with respect to the quantum fields r, tj, x, X + ) form of the classical action 
(4). After some algebra we get the following self-adjoint bilinear form 

( H TT H rr j Pl TX H TX + \ 
tt _ H m H r]X H vx + , 

11 IT IT IT IT ' \ V ) 



/ H TT 


H T rq 




H TX+ 


\ 


Hrq T 






TT 

vx + 




H X+r 


TT 

rL x + v 


TT 

rL x + x 


H x+x+ 


/ 


V H XT 




Hxx 


TT 

xx + 



with the following derivative structure of the operator H: 

/ D 2 + 2V^V^V V + + U QiD + + Q 3 \ 

The bilinear form of the action is matrix differential operator with fourth derivatives in the H TT 
sector and with second derivatives in the sector of usual scalar fields and in the mixed pieces. This 
structure of the bilinear form is similar to the one which is known from the theory of multiscalar 
GUT coupled to higher derivative gravity [26]. This is indeed natural, because what we are doing 
now is nothing but the study of induced quantum gravity (2) unified with the same GUT model. 
This analogy facilitates the calculations considerably, because of the following three reasons: 

i) The general expression for the divergent part of | Tr In H (f 0) is known from [26] 



,n — 4 

div 



ii n - 4 r r i „ . i „ „ i . . . . . 1 - - 2 

L 2 J d n x^ij tr | -P£Q 2 , + -P t Q 3 - -V/tPiQi ~ DP 1 Q 1 + -(P 1 Q 1 ) 



n 



+ h ( R ^ R ^ ~ I r2 ) ~ j^ R2 } + \ Tr ln { in + ^ AVa + t) } d%v + ( surface terms ) • 

The above formula includes the standard contribution of the second order operator 

^Trln {lU + E x V x + b\ 

2 ^ J div 

(see, for example, [f]), the contribution of the minimal fourth order operator first calculated in [fO] 
and also the contributions from the mixed sector derived in Ref. [26] for the case of higher derivative 
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gravity coupled to matter fields. The formula (11) shows, in particular, that we do not need explicit 
expression for N x . 

ii) Calculating the divergences with the use of the formula (11) one can learn the form of the beta- 
functions. However, since the formula (11) is the same as for the higher derivative gravity, it is easy 
to see that the general structure of the renormalization group equations in our theory is also the same 
as the one established in [26]. In particular, the quantum gravity corrections to the beta-functions of 
the nonminimal parameters £j are completely universal, they do not depend on the model and have 
the same form as for the one-scalar model [24]. Therefore all that we need is the form of the quantum 
gravity corrections to the beta-functions of the constants These corrections can be calculated 
on flat background and hence in what follows we put g^ v = T]^ v and switch off all the curvature 
dependent terms. 

iii) The structure of the expression (11) is direct generalization of the one we have already studied 
in [24] for the 1-scalar model. One can easily check that here, just as in the 1-scalar model, the 
contributions of trPiQi, tr (P1Q2 + P2Q1), tr (P1Q3 + P3Q1), and tr P2Q2 cancel, and therefore 
we need only the expressions for V^ v , U, P\, Qi, E X , D. Disregarding the curvature dependent 
terms one can obtain the following expressions for V, U, Q\, P\, D : 

V*" = -a 2 (Ci $ a $ a + C2 ¥>,V) 9*" , U = 0, 



Pi 



-101 



( 6$ a \ 



V 



D = H 



/-dZ5 ab + V ab 



v a 

3 



-C2Z8) + V/ V* 



(12) 



where Z = a 2 (Vcr) 2 + a(Dcr) and 



V ab = ^&5 ab + /i$ a $ 6 + h {^D rah D rcd + D rac D rhd \ <Z> c <Z> d + f 3 5 ab <p 2 + f 4 D rab <p+ f ^ Y <p* , 



X 



Vt = h$ a <p % + UD rac &\—)<p" 



Substituting these expressions into (11) we arrive at the explicit form of the corrections to the 
matter sector of T d i v from the quantum conformal factor 3 . 



,,n — 4 f ( \ "\ 



3 The renormalization in the vacuum sector was discussed in [24] for the case of the one-scalar model. Its form 
doesn't depend on the number of scalar fields and that is why we do not discuss it here. 
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+2« 4 GC 2 (1 - Ci " C2 + CiC2)$Vv? + lf2<* 2 <;H* a D abc <f> b ) 2 + ha 2 [C 2 + 4CiC 2 + c 2 2 ] $V + ¥> 

+Ua 2 [C2+4C1C2 + C2 2 ] ($ a D abc $ b )<pt (y) V' + 6/ 5 « 2 C 2 2 ^ + ^ 2 | , (13) 

where e = (47r) 2 (ra — 4) is the parameter of dimensional regularization, and we disregarded all surface 
and matter-independent terms. 

One can see that there are no any divergences which lead to renormalization of scalar fields $, (p. 
Thus the only modifications due to the contributions of the conformal factor are in the renormaliza- 
tion of couplings /i,2,..,5 and parameters £^2, and the renormalization group equations include only 
(3— functions, but not 7-functions. The contribution of the conformal factor to the effective potential 
of the scalar fields depends on the /3-functions only. 



4 Renormalization group equations 

Now we are in a position to perform the renormalization group study of the back-reaction of vacuum 
to the matter fields. For our purposes it is convenient to use the formulation of the renormalization 
group in curved space-time, given in [17, 1]. The renormalization group equation for the effective 
action in curved space-time has the form 

{^ + ^T P -J d ^ 91w s^xj } r[w ' p ' 9 ^ » ] = ' (14) 

where w is the full set of the quantum fields (gauge, spinor, Higgs and conformal factor a) and 
p - complete set of parameters including £'s. The solution corresponding to the desirable scaling 
behavior has the form 

r[e" 2 *^, w, p, fi] = Tlg^, w(t), p(t), fi] , (15) 

where fj, is the dimensional parameter of renormalization. Effective fields and coupling constants 
obey the equations 

(47r) 2 ^ = ( 7u , + ^)^, w(0) = w, (16) 



(An) 2 ^ = (3 p +pd pi P(0)=P, (17) 

where 7 and (3 functions are defined as usual. According to the results of the previous section the 
renormalization of the fields and all couplings except fi and £j are not modified by the contributions 
of the quantum field a. Thus the /3-functions for the gauge and Yukawa coupling constants are just 
the ones derived in [25]. On the other hand, these contributions to the fij t and (3^ all have universal 
form and do not depend on the gauge group of the theory. Thus we find for our theory (from this 
moment we will be using the variables Q = 1 — 6£j for convenience) : 

P h =pf + W h , ft-^ + Aft., (i hk =(i [ S, (i g = (if\ (18) 

where f3^ is the /3-function for the effective parameter p in curved space-time without back reaction 
of vacuum (or other form of quantum gravity) and A(3 p are the quantum gravity corrections. In our 
case, contrary to the high derivative gravity [26] those A/3 p 's are nonzero only for the scalar and the 
nonminimal parameters. 
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According to recent communications (see, for example, [28]) the values of £j are very important 
at the energies between the Planck scale Mpi = W 19 GeV and the unification scale Mj = W 14 GeV, 
because there is a hope to meet natural inflation for this rate of character energies. Thus our purpose 
is to evaluate the running of £j (or Q) backward from the high energy Planck scale. As a lower end of 
the energy interval one can take the unification point Mj = W 14 GeV but, since it doesn't lead to the 
serious changes in the calculations we shall take, as a lower limit, the Fermi scale Mp = 100 GeV. 

4.1 SU(N) model: running away from the finite theory 

As an examples of the effect of the quantum conformal factor we shall consider two toy models 
with two scalar fields: one which possesses the one-loop finiteness without supersymmetry [25], and 
another one which admits N = 2 supersymmetry but is not finite. 

Let us start with the generic model (4). The counterterms (13) lead to the following expressions 
for A/V 

(4K) 2 Af3 fl = 12« 2 /i Ci 2 + 4« 4 Ci 2 (Ci - I) 2 , /i(0) = fi ; (19) 



(4tt) 2 A/3 /2 = 12« 2 /2Ci 2 



/ 2 (0) = / 2 ; (20) 



(4tt) 2 A/3 /3 = 2« 2 h (Ci 2 + 46C2 + C 2 2 ) + 4« 4 CiC 2 (1 - Ci - C2 + C1C2) , / 3 (o) = h ; (21) 



(4tt) 2 A/3 /4 =2« 2 / 4 (Ci 2 + 4CiC2 + C: 



2) 1 



(4tt) 2 A/3 /5 = 12a 2 h Q + 4« 4 C 2 2 (C 2 - I)' , 



/ 4 (0) = / 4 ; (22) 
/ 5 (0) = / 5 (23) 



and 



(4tt) 2 Aft . = 2« 2 C| (0 - 1) , C'(0) = Ci • (24) 

Let us now consider the asymptotic behavior of the effective couplings fi(t),Q(t). As it was 
already mentioned above, all the /3-functions (corresponding to gauge, Yukawa, scalar couplings and 
to the nonminimal parameters £) vanish in the conformal fixed point 

h k (t) = h* k , Mt) = f*, = 0, (25) 

where /*, h\ are the values corresponding to the fixed point in flat space-time. The corrections 
A/3 indicate that there is also a second "minimal" fixed point with the same (25) solutions for h k {t) 
and fi (t) but with Q = 1 (this corresponds to the £j = in (4) , that is why we call this fixed point 
minimal) . The behavior of the effective charges and effective potential in the vicinity of the minimal 
fixed point has been studied in [24] for the one-scalar model. Now we are interested in the behavior 
of the effective charges close to the conformal fixed point. 

As far as the renormalization group equations for g(t), hk(t), fi(t), £?'(£) are ver y cumbersome, 
it is reasonable to take particular models in which the study of these equations performs easier. 
As a first example, let us take a particular model of (4) with N = 8 and one scalar multiplet in 
the adjoint representation. When the number of spinor multiplets is m = 84 the theory is one-loop 
finite in the flat space-time (see Ref. [25] for a detailed discussion of such a model without taking 
the quantum effects of the conformal factor into account) . To perform the numerical analysis of the 
renormalization group equations one has to impose the initial conditions at far UV. We suppose that 
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the running from the conformal fixed point happens because of the two-loop contributions which 
violate the £ = g condition and thus switch on the interaction with the conformal factor. Thus, to 
impose the initial conditions one has to evaluate the two-loop effects. The two-loop contributions 

4 

are, generally speaking, proportional to j^pp multiplied by the combinatorial factor. Let us take the 
typical value of the parameter ( as 10~ 3 , while the coupling constant is g 2 = O.f . Then the numerical 
analysis enables one to trace the behavior of the effective couplings ((t), /i,2 (i) for 
> t = ln(fi/Mpi) > —39 that corresponds to the running between fi = Mp\ = W 19 GeV and 
fi = M F = W 2 GeV. 

We have performed the numerical analysis of the full system of the renormalization group equa- 
tions for the effective gauge, Yukawa, scalar couplings and ((t) with taking into account the cor- 
rections (19) - (24). As an initial point of the renormalization group flow in far UV we choose the 
values (see [25] for the details) 

h 2 = ; ft « 0.43621 , / 2 * « 0.30277 

that provides the one-loop finiteness in the theory without back reaction of vacuum. The results for 
the contribution to the couplings due to back-reaction, that is the differences Dfj and D(j between 
the couplings in the with and without the back-reaction of vacuum, are plotted in Fig. 1. It is easy 
to see that the plots for the deviations look like a plane lines. The reason for this is that the effect of 
the back reaction is very week and in such a "small" interval on the logarithmic scale the non-linear 
function looks like a linear. Also the numerical values presented at the Fig. 1 show that the effect is 
very small. For instance, the deviation of fi^ from the finite fixed point is about six orders smaller 
than the value of the coupling itself in this fixed point. One has to notice that since the effect is 
almost linear, it is quite easy to construct the same plots for other initial deviations DQ\. 

The similar numerical analysis has been performed for the more interesting SU(5) model (4) - 
(6) with two kinds of scalar fields and with 

m=l, ra=15; ^1,2,3^) = ^1,2,3 , fi,...,5 (*) = /i,...,5 I 

where 

h\ 2 3 = (1.421 , 1.681, 2.361) ; /*,2,3,4,5 = (0-659, 1.293, 0.324, 1.677, 1.039). 

The results of the numerical solution are presented at the Fig. 2 .The parameter D(i(0) here is taken 
0.01 as in the previous example, and we kept the constraint D(i = DQ2 for the sake of simplicity. 
Qualitatively the results are the same as in the one-scalar case - the effect exists but it is tiny. 




Figure 1: Contribution to the running coupling constants due to back-reaction. SVfS] model with one scalar 
field and 84 spinor multiplets. 
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Figure 2: Contribution to the running coupling constants due to back-reaction. SV^] model with two scalar 
fields and with the spinor content m = 1, n = 15. 



4.2 Model with the rigid breaking of supersymmetry 

Let us now consider the SU(2) gauge model with an action [29, 16, 1]: 




+i if) a 7^ V ab tp b - hitpa F abc (p b tp c -i h 2 il> a j5 F abc A b lj) c + 



-\htf) 2 -\m 2 f -\h<?*? -\h [p 2 ^ 2 - (V^A) 2 ] } 



(26) 



where ip and A are scalar and pseudoscalar fields in the adjont representation of the SU(2) gauge 
group. Other notations are obvious. In flat space, for the particular values of the couplings 



the model (26) possesses the N = 2 supersymmetry. For other values of the couplings one meets 
the renormalizable theory which has also second non-supersymmetric renormalization group fixed 
point [29]. In curved space time this model is also renormalizable (if only the nonminimal terms and 
the action of vacuum are introduced), and admits the asymptotic conformal invariance for (27) and 



Taking into account our previous results we arrive at the renormalization group equations for the 
above model with the back reaction of vacuum. The gauge coupling behaves as [29] 



and supersymmetry fixes Yukawa constants to be h\ = h\ = g 2 . In flat space-time those relations 
hold under renormalization [29]. The same happens when we put our theory in curved space-time and 
take the back reaction of vacuum into account, because the quantum conformal factor contributes 
to the /3-functions of the scalar couplings and nonminimal parameters. We are interested in the 
renormalization group flow from the supersymmetric and conformal invariant UV fixed point, that 
is why we have to impose the constraints (27) on the initial data. As far as the back reaction doesn't 
concern the behavior of the Yukawa couplings, the relations (27) for the Yukawa couplings hold for 
all scales and they can be indeed used in the RG equations directly, while /i,2,3,4 should be regarded 



fi = f2 = h = 0, U = hi = h 2 2 = g 2 



(27) 



6+6 = 1/6. 
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as arbitrary quantities. After we use the relations between gauge and Yukawa coupling constants, 
the RG equations for /i,2,3,4(i) become 

(47T) 2 ^ = 11/2 + 3/ 2 + 8/s/4 + g/ 2 _ 8fig 2 _ 8g 4 + l2 0? hfi + 4^ (Cl " l) 2 , 

(4tt) 2 ^ = 11/| + 3/| + 8/3/4 + 8/ 4 2 - 8f 2 g 2 - 8g 4 + 12« 2 / 2 C 2 2 + 4« 4 C 2 2 (C 2 - I) 2 , 

(4vr) 2 ^ = 4/| + (5/3 + 4/4) (A + / 2 ) + 8/ 4 2 - 8/3/ " 8g 4 + 2« 2 / 3 (d 2 + K1C2 + C 2 ) + 
+ 8a 4 CiC 2 (l-Ci - C2 + C1C2), 

(4tt) 2 ^ = 6/ 2 - 8/4/ " 6</ 4 + 8/3/4 + 2/4 (A + / 2 ) + 2« 2 / 4 (d 2 + 4dC 2 + Cl) , (28) 



while for £1,2 (i) we meet the equations: 



(4vr) 2 ^ = Ci (5/i - 4 5 2 ) + C 2 (3/ 3 + 4/4) + 2« 2 C 2 (Ci - 1) 



(47T) 



dt 



C2 (5/2 - 4 5 2 ) + Ci (3/s + 4/4) + 2« 2 C 2 2 (C 2 - 1 



(29) 



In order to study the renormalization group behaviour of the effective couplings f(t), £(t) we 
have to choose the initial value of g 2 and the initial variation _D(j = — DQ2 / 0. As in the previous 
case we take in the far UV g 2 = 0.1 and (j = 0.001 . The results of the numerical analysis, 
concerning the contributions due to back-reaction, are plotted in Figs. 3 for the case (j = 0.001. 





Figure 3: Contribution of back-reaction to the running couplings fj(t) and Cj(t)- First example with 

d(0) = -c 2 (o) =0.001. 



One can see that in this case the effect of quantum conformal factor is essentially the same as in 
the previous case of the SU(8) and SU(5) models. Because of those effects one meets the nonzero 
/i,2,3 couplings with the first two having the values about 10~ 8 at the Mj unification scale and 
about 10~ 6 at the Mp Fermi scale, while they remain identically zero without the back reaction of 
the conformal factor. 

In order to give some speculation on the strongly coupled theories we also investigated the case 
with initial values of g 2 = 0.5 and £ = 0.1. The results are plotted at Fig. 4. 
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Figure 4: Contribution of back-reaction to the running couplings fj(t) and CjW- 9 2 = 0-5 and Ci(0) = 

-C 2 (o) = o.i. 




Figure 5: Contribution of back-reaction to the running couplings fj(t) and Cj(t)- 9 2 = 0-5 and Ci(0) = 
— £2(0) = 0.5. This example was included only to illustrate the nonrealistic order of magnitude for g 2 and 
initial deviations DQ for which the nonlinear behavior shows up. 

5 Conclusion 

The arguments were presented that in four-dimensional quantum field theory in curved space-time the 
local conformal symmetry can not be exact, and that it can be only approximate symmetry which 
serves as the initial point for the renormalization group equations in the far UV limit. Starting 
from the conformal invariant theory one meets the trace anomaly and consequent propagation of the 
conformal factor of the metric. We have considered the interaction between the quantized conformal 
factor and the matter fields in a theory with approximate conformal symmetry, in a region close to 
the scale of asymptotic freedom and asymptotic conformal invariance. 

The contributions of the conformal factor to the /3-functions of the scalar coupling constants have 
been calculated and it was shown that these contributions drive the theory out from the conformal 
point and also modify a values of the scalar coupling constants. If the starting model in far UV 
possesses supersymmetry and (or) finiteness, the interaction with the conformal factor leads to the 
violation of those properties at lower energies. For the initial violations about DQ = 0.001 the 
numerical values of the deviations of the scalar couplings from the symmetric state range from 10~ 8 
to 10~ 6 at the Fermi scale and about one order less at the Unification scale. Indeed the effect is 
very weak and therefore all the dependencies are very close to be linear. In particular the above 
deviations may be considered as proportional to the squares of the initial violations (DQ) 2 of the 
conformal invariance. At higher loops the supersymmetry breaking in the scalar sector should violate 
supersymmetry in the Yukawa and gauge interactions. As a result the relations between the masses 
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of the particles at low energies may differ from the ones which are expected in supersymmetric 
GUT's and the supersymmetry is not seen at low energies. However, since the numerical effect of 
such a violation will be extremely small, it is very difficult to see how one can observe it. One can 
say that if the local conformal invariance exists as a high energy symmetry, it holds as a very good 
approximation at lower energies. 
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